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Abstract 

A mixing chamber used in rocket engine combustion testing at NASA Stennis Space Center 
is modeled by a second order nonlinear MIMO system. The mixer is used to condition the 
thermodynamic properties of cryogenic liquid propellant by controlled injection of the same 
substance in the gaseous phase. The three inputs of the mixer are the positions of the valves 
regulating the liquid and gas flows at the inlets, and the position of the exit valve regulating the 
flow of conditioned propellant. The outputs to be tracked and/or regulated are mixer internal 
pressure, exit mass flow, and exit temperature. The outputs must conform to test specifications 
dictated by the type of rocket engine or component being tested downstream of the mixer. 
Feedback linearization is used to achieve tracking and regulation of the outputs. It is shown 
that the system is minimum-phase provided certain conditions on the parameters are satisfied. 
The conditions are shown to have physical interpretation. 


1 Introduction 

The NASA John C. Stennis Space Center (SSC) conducts extensive ground-based testing and flight 
certification of rocket engines, in particular, of the Space Shuttle Main Engine (SSME). Combustion 
chambers and turbomachinery related to rocket engines are also tested at SSC. This work is part of an 
on-going effort to develop a software package providing flexibility in simulation and control tasks [10, 
2, 3] frequently found in test operations at SSC. Test conditions require that liquid propellants, 
namely liquid oxigen and liquid hydrogen (LH2) be supplied to the engine or component at very 
precise conditions of temperature, pressure and mass flow rate. An excess or deficiency in any of 
these three flow parameters may result in damaged components or in a sub-optimal test. To achieve 
the required conditions, the delivery system includes a mixing chamber, henceforth referred to as 
“mixer” , The mixer subsystem is depicted in Fig. 1. LH2 is stored in the run tank, which is kept at a 
constant pressure by an independent control loop. Gaseous hydrogen (GH2) is stored in high pressure 
bottles, and has a higher temperature than the LH2. One control valve, referred to as the “liquid 
valve” is used to manipulate the flow of LH2 into the mixer. There is also a “gas valve” and an “exit 
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Figure 1: Mixer Subsystem 

valve” . The positions of the valves constitute the only control variables for the mixer subsystem. 
The objective of the mixer control system is to specify a MIMO controller which achieves tracking 
and regulation of mixer outputs to desired values. The outputs of interest are the mixer pressure, 
exit flow temperature, and exit mass flow. The MIMO system, thus, has three inputs and three 
outputs. A dynamic model and controller are intended to replace the current method of operation 
of the mixer, which uses only a steady-state thermodynamic model and extensive heuristics. The 
use of the modeling and control techniques described herein is expected to provide great flexibility 
and better mixer performance. 


2 Transient Thermodynamics Model of the Mixer 

2.1 Valve Models 

The flow wi of LH2 through the liquid valve depends on the pressure difference across the valve, the 
density of the source LH2 and the liquid valve opening coefficient C v i according to [1]: 

w t = Ci<W(fi - P)Pi (!) 


In English units, Pi is the source LH2 pressure and P is the mixer pressure in psia; pi is the source 
LH2 density in lbm/ft 3 , and C\ = 1.76 x 10 -2 . For gas flow, there is a separate expression [1]: 


w g = 



,when P g < 2 P 


( 2 ) 


\ C±C V g^fTgp 


,when P g > 2 P (choked flow) 
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where T g is the source GH2 temperature in degrees Rankine (°i?), P g is the source GH2 pressure 
in psia, and C vg is the gas valve opening coefficient. In English units, C 2 — 2.857 x 10“ 2 and 
C 4 = 2.423 x 10" 2 . 


2.2 Mass and Energy Balances 

The set of two nonlinear differential equations which constitute the model can be derived from the 
conservation of mass for compressible flow and from the First Law of Thermodynamics for a control 
volume including transient terms [11]. The equations may be expressed in terms of internal energy 
and density derivatives as: 


p = ~(wi+w g -w e ) 

u = -j- [wi(hi - u) 4* Wg(h g “ u) - w e (h e - «)] 

r ' 


(3) 

(4) 


where V is the fixed mixer volume, hi, h g and h e are the source liquid, source gas, and exit enthalpies, 
respectively. The mass flows Wi are functions of pressure, and are calculated according to Eqs.( 1) 
and ( 2). 


2.3 State Variable Model 

Further manipulations are required to put the above physical model in standard state-space form. 
It is assumed that enthalpy is conserved across the exit valve, and the definition [11] fr = u + CP/p, 
where C is a constant depending on the choice of units, is used. For English units, C = 0.185. The 
resulting model is compactly expressed as 

P = fi(P(p,u),p)C v i + f 2 (P(p,u),p)C vg + f 3 (P(p,u),p)C ve (5) 

u = g 1 (P(p,u),p,u)C v i + g 2 (P(p,u),p,u)C vg + g 3 (P(p,u), p,u)C ve (6) 


where 


fi 

h 

h 

f 3 

g i 

P2 

gz 


cW{p 1 -p)pi 

C' 2 , if Pg > 2 P 

C -4 \J Pg ~ P 2 , if Pg < 2P 

hi - u' 


f l 

h 

c'i 


hn 


lP-PsP 


The constants are given by C[ = Ci/V, C' 2 = C 2 \fP) j Pg/V , C' z = -Cz/V = -Ci/V , C z = CCi 
and C'i = Ci^/T g p g /V. 
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2.4 Output Definitions and Controls Model 

Mixer operation requires the simultaneous tracking of exit temperature, exit flow, and mixer pressure. 
Exit temperature and mixer pressure are functions of the energy and density states. Computation 
of these functions requires the intervention of thermodynamic tables and interpolation algorithms 
which cannot be represented in closed form. However [7], to each exit temperature and mixer 
pressure combination in the expected range of operation there corresponds a unique value of the 
state [p, u] T . Therefore, it is convenient to specify density and energy as outputs, along with 
exit mass flow. Desired state trajectories that correspond to desired exit temperature and mixer 
pressure can be either pre-computed off-line, or they can be obtained through a reference prefilter. 
With these output definitions, the controls model becomes 


< 


p 

it 


fi(p,u) /2O, u) 

9i(p,u) 9i{p,u) 


U = Ip, u, -v hC ve ] T 


h{p, u) 
9s{p, u) 


1 

Cvl 


C V Q 


Cye 


where [C v i C vg C ve ] T is the control vector and y is the output vector. 


(7) 


3 Model Properties 


3.1 Uniqueness of Equilibrium 


The equilibrium point indicates the steady values of the density and internal energy of the mixer 
when the control valves are set at fixed positions, for a given set of input flow characteristics. 
For given values of the input fluid properties and C v coefficients, setting p cv = 0 establishes that 
rh e = trig + mi and results in an expression relating the equilibrium density p cv to the equilibrium 
mixer pressure P : 

Pcv = p{P ) ( 8 ) 


Setting u cv = 0 gives the equilibrium exit enthalpy h e = h cv in terms of the input enthalpies and 
the mass flows: 

_ - h ( P ) ( 9) 


hcv 


rhihi + rhghg 


mi + m g 

The above enthalpy must match the thermodynamic property data at P and p cv of Eq.( 8), that is 


hcv — hth ( Pcv ) P ) 


Substituting Eqs.( 8) and ( 9) into the above equation results in a single expression which gives the 
equilibrium pressure: 

h(P) = h th (p(P),P) (10) 

A graphical interpretation of the equilibrium solution is shown in Fig. 2. The curve pP vs. P has 
been drawn on the base plane. This plane curve is mapped to a space curve by the thermodynamic 
property function hth • 'The equilibrium exit enthalpy h(P) is a function of pressure only and the 
corresponding surface has also been graphed. The point where the space curve pierces the surface is 


tu/lbm 


2003 American Control Conference , Denver , June 4 - b. 


2500 


2000-4 


1500-4 


Graphical interpretation of Mixer Equilibrium 


Steady-state enthalpy 
surface 


( Enthalpy from 
thermodynamic 
data 





2GG3 American Control Conference, Denver, June 4-6 . 


6 


the equilibrium point of the system. The monotonicity of the curve and surface indicate that there 
exists only one equilibrium point in the range of interest of actual mixer operation. An iterative 
procedure is developed in [7] that calculates the equilibrium point, given thermodynamic parameters 
and fixed control inputs. 


3.2 Valve Positions for a Prescribed Operating Condition (Reachability) 

Having three independent controls in a two-state model allows the selection of steady values for the 
states, and, in addition, an extra degree of freedom is available. This degree of freedom can be used 
to fix the exit mass flow with the desired thermodynamic properties, as shown next. Supppose it is 
desired to have a given exit mass flow, with prescribed temperature (measured at the outlet of the 
exit valve). Let these quantities be denoted by w e and T e . The back pressure P s at the outlet valve 
exit is assumed to be constant. The problem is to determine the valve coefficients that achieve this. 
Two degrees of freedom are used for w e and T e , and the third one is used to meet a desired mixer 
operating pressure, P. The back pressure P s and T e determine the enthalpy h, which is constant 
across the valve. Therefore, the enthalpy and mixer pressure at the exit valve inlet are known 
and, in turn, determine the mixer density from thermodynamic data. A static energy and material 
balance gives the required input flows that achieve a prescribed flow-exit enthalpy combination. The 
following formulas are straightforward to derive: 


wi 


w 


9 


W e (h — hg) 

h[ hg 

w e (hi - h) 
hi hg 


( 11 ) 


where h g and hi are the gas and liquid supply enthalpies, respectively. Using the exit flow and 
the input equilibrium flows from Eq.( 11), the three valve coefficients are determined from Eq.( 1) 
and ( 2). This procedure is complementary to the one described in the previous section, which 
calculates the equilibrium point. It shows, in a rather informal way, that the system is reachable in 
the thermodynamic range of interest. 


4 Feedback Linearization Analysis 

Feedback linearization [12, 8, 5], has been succesfully applied in a variety of cases related to thermal 
process control; for instance in [6, 9, 4]. In this section the technique is applied to the mixer. The 
model of Eq.( 7) has a characteristic which prevents direct application of input/output linearization 
theory. It is seen that the third component of y is statically related to one of the control inputs, 
namely C V€ . This implies that the partial relative degree of y% is zero. Conceivably, one could use 
a time-explicit control of the form 

C ve(t,p,u) — y h{ p Mp) 

to attain perfect tracking of the output mass flow rate and use the two remaining controls to force 
the flow to have the desired pressure and temperature. This approach, however, is not robust. 
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Therefore a a feedforward law is not considered. One way to get around this problem altogether is 
to augment the exit valve channel with an integrator, that is, let: 

C ve = v 


where v is a new control input. Now C ve is regarded as a state, and the resulting system is of third 
order, with three inputs and outputs. If the arguments of functions gi and fi are dropped from the 
notation, the new system equations become: 


p — flCvl + f 2 C V g + fs C ve 

(12) 

u = 9iC v i 4* 92C vg -f gzC ve 

(13) 

Cye — V 

(14) 

y = [p u -Vf 3 C ve ] T 

(15) 


Upon differentiating the outputs once and rearranging, the derivatives can be written as: 

y = D + Ew 

where y = [y x y 2 y^} 7 , w - [C„i C vg v] T and 


D = 


fs^ve 

9%Eve 



E = 


fi 

9i 


h 

92 


-VC ve [^(^fi + f^ 5 i) + a^/i] -VC ve [^(^f 2 + §£52) + 7^/2] 


0 

0 

-Vh 


It is seen that the partial relative degrees are all one. Provided E is invertible in a region fi of the 
state space, the feedback law 

w = E- 1 (y d -T(y-y d )-D) (16) 


achieves exact linearization of the system, with tracking error dynamics given by 


(y - yd) + r (y -yd) = 0 


If T is chosen as a diagonal positive-definite matrix, the resulting control law is a “decoupling 
control” , since the dynamics of the tracking errors are noninteracting [5]. If the appropriate function 
definitions are substituted, the form of the E matrix is as follows: 


E = 


c{ s/iP, - P)PI 
C 1 \A P J - 


VC ve C'. 




2 s /(P~P s )p 




v C-ue /2 ^3 
2s/(P-P s )P 


f2 



[off + 


0 

0 

-VC' S ^(P - Ps)p 


4.1 Input-Output Linearizability of Augmented Model 

The ability to construct a feedback linearization controller hinges on stable zero dynamics of the 
augmented system, and on the invertibility of matrix E. 
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4.1.1 Invertibility of E 

By inspection, it is readily seen that the first two rows of E are linearly independent provided h g ^ hi. 
This has a direct physical interpretation: if the two fluids have the same thermal properties (i.e., 
enthalpies), the ability to change the thermal properties of the mixture by changing the relative 
flows is lost. Gas and liquid enthalpies are different for the expected operating conditions. The 
third row is linearly independent from the first two if P — P s >0, which is also true for the mixer. 
Therefore E is invertible over the whole range of expected mixer operating conditions. 


4.1.2 Zero Dynamics of the Augmented Model 


As it is known [5], the zero dynamics is unchanged upon augmentation with an integrator. This 
implies that one may examine the non-augmented model for zero dynamics and draw conclusions 
about the augmented model’s zero dynamics. Suppose the exit flow is to be held constant at a value 
I 30 . The only way in which this can be achieved is by letting 

- - v m> (17) 

at all times. Differentiating the other two outputs and equating them to zero results in 

C v i(t) = (-h{t)C vg {t) + ^) (18) 

c ” (t) = “So( 9,wc *' (!) + w) (19) 

Upon substitution and rearrangement, it is seen that there exist three uniquely defined control inputs 

which hold the outputs constant. Since two of the outputs coincide with the states, it trivially follows 

that they are kept bounded, and therefore the system has stable zero dynamics (i.e., the system is 
minimum-phase). The inputs are given by Eq. 17 and, dropping the time notation: 


C vg — 
Cvi — 


Xmh. ( 23 „ 

y y/3 h) 

fi92 - 9ih 


Yz o_ 
Vfi 


1 - 


/ 1/2 f 9s _ £1 
fi 92 - £ 1/2 V /3 h 


( 20 ) 

( 21 ) 


Note that the above formulas can be used to find the valve positions at which the system has a 
prescribed outflow and a pair of thermodynamic properties. Therefore it produces the same results 
as in Section 3.2. 


5 Simulation Study of the Feedback Linearization Controller 

Functions that compute the pressure and its partial derivatives at each time step were implemented 
in Matlab. Table 1 summarizes the parameters used in the simulations. The decoupled control 
gains were set as 10 for the density, 10 for the energy and 5 for the exit flow. The controller has 
exact knowledge of all parameters, except the initial conditions. The objective in Simulation 1 is 
to transfer the mixer from the initial conditions to a pressure of 6000 psia, exit temperature of 
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Parameter 

Value 

Units 

Pi 

8500 

psia 

Ti 

-340 

°F 

P 9 

13500 

psia 

T g 

90 

op 

Ps 

5533 

psia 

Po 

3.85 

lbm/ffc 3 

U 0 

180 

Btu/lbm 


Table 1: Feedback Linearization Parameters 

-200°F and exit flow of 40 lbm/s. With the known and fixed back pressure it is determined that 
the mixer states must be transferred to p — 2.93 Ibm/ft 6 and u = 400.9 Btu/lbm for the pressure 
and exit temperature to converge to their desired values. The reference states and flow are specified 
as constants, leaving the trajectories to be determined by the controller. In Simulation 2, the mixer 
is asked to reach the previous equilibrium point and then raise the exit temperature to 0°^ and 
the pressure to 7000 psia, while reducing the exit flow to 10 lbm/s. This corresponds to lowering 
the density to p = 2.143 and raising the energy to u = 995.9. This time, the states and the flow 
are required to track ramp functions which transfer from old to new values in 5 seconds. Figs. 3 
and 4 show the results of Simulations 1 and 2, respectively. The steady values required are exactly 
achieved. However, the transient response of the exit flow in the simulations has some overshoot. 
This may be reduced by modifying the control gains, but perhaps at the expense of slowing down 
the response. 

6 Conclusions 

A mathematical model for the mixer is presented. The model consists of a second order nonlinear 
MIMO system having density and internal energy as states. The inputs are valve positions repre- 
sented by flow coefficients, and the controlled outputs are mixer pressure, exit temperature and exit 
mass flow. Model reliability is ensured by the use of function calls to real thermodynamic data. It is 
shown that, in the expected range of operation, the mixer has a single equilibrium point for each set 
of fixed valve coefficients. Conversely, the reachability of operating conditions is shown by the pos- 
sibility of determining a unique set of valve coefficients that achieve a desired set ot steady outputs. 
A feedback linearization controller is derived for the mixer system, and it is shown that the system 
is minimum-phase and the linearizing matrix is invertible under parametric conditions having nice 
physical interpretations. The feedback linearization controller gives good results. Further work in 
the topic includes examining other control techniques providing robustness against variations in the 
supply fluid properties (bottle and tank depletion) when they cannot be measured. 
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